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Abstract 

We introduce and discuss a new class of solutions of the DMPK equation 
in which some of the eigenvalues are grouped into clusters which are conserved 
in the asymptotic large distance limit (i.e. as the length of the wire increases). 
We give an explicit expression for the asymptotic expansion of these solutions 
and suggest some possible applications. In particular these new solution could 
be useful to avoid the quasi one dimensional constraint in the DMPK equation 
and to study the crossover between the metallic and insulating phases. 
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1 Introduction 



One of the most interesting tools to describe the electron transport properties of 
quantum wires is the so called Dorokhov Mello Pereyra Kumar (DMPK) equa- 
tion [1]. This equation describes the evolution of the joint probability distribution 
of the transmission eigenvalues P({Aj}) as the length of the wire increases and has 
been the subject of intense study in these last years [2]- [11]. Among the several 
remarkable features of this equation, the most interesting one is that it can be 
exactly mapped onto the radial part of the Laplace-Beltrami operator of suitable 
symmetric spaces [2]. Thanks to this mapping it is possible to write the Green 
function of the DMPK equation in terms of the so called Zonal Spherical Functions 
(ZSF) of the corresponding symmetric spaces. In the [3 = 2 case the ZSF can be 
written explicitly in terms of ordinary ipergeometric functions [12] thus leading to 
an exact solution of the DMPK equation [4] while for /3 = 1,4 one can rely on a 
powerful asymptotic expansion due to Harish-Chandra [13]. This approach leads 
to a rather involved expression for P{{Xi}) which however drastically simplifies in 
the two limits of very short wires (metallic regime) and very long wires (insulating 
regime), thus allowing to evaluate all the quantities of interest (see [4] and [5] for a 
detailed discussion). 

However, despite these remarkable results there are two major drawbacks in the 
DMPK approach to quantum wires. The first is that the DMPK description only 
holds in the (quasi) one dimensional limit. The second is that the above mentioned 
solution (even in the simpler [3 = 2 case) does not allow to study the intermediate 
region between the metallic and the insulating regimes (the so called "cross-over 
regime") where no simplifying approximation is allowed. This region has recently 
attracted much interest since both numerical and analytical results seem to indicate 
a rather non-trivial (with possibly a non-analytic point [6]) behaviour of P({Aj}) 
exactly in this regime [7, 8, 9]. 

In these last years several interesting approach have been proposed to over- 
come these problems. Various generahzations of the DMPK equation have been 
suggested [11] [10] to avoid the quasi one dimensional limit. However in all these 
generalized equations most of the nice properties of the DMPK equation are lost 
(mainly due to the fact that the description in terms of symmetric spaces is lost) 
and only few informations on the expected joint probability density of the transmis- 
sion eigenvalues can be obtained. The aim of this paper is to propose a completely 
different approach to these problems, which instead fully exploits the power of the 
symmetric space structure which is behind the DMPK equation. Instead of mod- 
ifying the DMPK equation , we shall keep it unchanged, but shall look for a set 
of special solutions (with non trivial initial conditions) which break the isotropy 
ansatz. To this end we shall use the fact that the DMPK equation can be mapped 
into the evolution operator of a class of ID quantum integrable models known as 
Calogero-Sutherland (C-S) models (for a review see [14]). As a consequence of their 
exact integrability it is possible to show that in these models, besides the well known 
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symmetric solution, a wide class of non trivial (but exact) solutions exist in which 
the particles are grouped into clusters which survive in the asymptotic limit. The 
clusters of particles become, once the mapping with the DMPK equation is per- 
formed, clusters of eigenvalues. The exact integrability of the C-S model ensures 
that this asymmetric distribution of eigenvalues survives in the asymptotic limit 
and the remarkable properties of the underlying symmetric space allow to write 
explicitly such asymptotic expansion. 

This paper is organized as follows: In sect 2 we shall recall some known results 
on quantum transport, focusing on the DMPK equation in order to fix notations. 
In sect 3 we shall review the main concepts associated to C-S models and the 
mathematical tools applied in order to find their solutions. In sect 4 we shall 
describe the new class of solutions of the DMPK equation and in sect 5 we shall 
give some hints on how these new solutions could be used in order to address the 
open problems mentioned in the introduction. 



2 Quantum transport in a wire and DMPK equa- 
tion 

A mesoscopic conductor can be modeled as a disordered region located between 
two ideal leads connected to two electron reservoirs; at very low temperatures, 
frequencies and voltages the scattering phenomena inside the wire are supposed 
to be elastic, and the electron wave functions are assumed to be phase coherent. 
Given a quantum wire of lenght L and width W {L ^ W) it is possible to find a 
finite but large number N of Fermi channels associated to the conduction electrons. 
Adopting the notation used in [1], we call / and O the A'"-dimensional vectors 
associated respectively to the amplitudes of the incoming and outgoing waves on 
the left of the wire and /' and O' the corrisponding vectors on the right. 

The scattering phenomena inside the wire can be described by means of the 
2iV-dimensional Transfer matrix M, which connects the amplitudes of each channel 
located on the left to the ones on the right of the wire: 

(1) 

The matrix M can be decomposed in a radial and an angular part according to the 
polar decomposition [1]: 

with Ui unitary matrices NxN, which assume the role of angular coordinates and 
A real diagonal matrix with non negative elements {Ai . . . Ajv}, as radial coordinates. 

In [1] the general symmetries of the system (time reversal and/or rotational 
symmetry) fix the matrix M to a particular Lie group: 
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• if time reversal symmetry (TRS) is present and there is no spin-rotational 
symmetry (SRS) inside the wire M e Sp {2N,R); 

• if time reversal symmetry and spin-rotational symmetry are both present in- 
side the wire M e SO*{4:N); 

• if time reversal symmetry is absent M G SU {N, N) . 

Experimentally these three symmetry classes correspond to the presence or absence 
of a magnetic field (which breaks the time reversal symmetry) and to the relevance 
of the spin-orbit interaction term (which breaks spin-rotational symmetry) inside 
the wire. 

The general symmetries of the Transfer matrix lead to a possible formulation 
of the problem by means of the tools offered by Random Matrix Theory; each 
symmetry class is assigned to an ensemble of tranfer matrices, associated to a 
parameter f3 which can be interpreted as the number of degrees of freedom of each 
matrix element. 

In the framework of the transfer matrix approach to qiiantum wires the DMPK 
equation [1] describes the evolution of the probability distribution Pl of the trans- 
mission cigcnparamctcrs {Aj} as a function of the reduced Icnght s = L/l {L being 
the length of the wire and / the mean free path associated to the electron): 

BP 2 ^ f) r) 

with 7 = /3A^ + 2-/3 

where is the number of Fermi channels associated to the conduction electrons, 
j3 = {1, 2, 4} is the parameter associated to the symmetry of the system, is the 
Jacobian from the space of the whole matrix M to the space of the eigenparameters 

TV N 

^/3({A.}) = n n iA.-A,i^ 

i=lj=i+l 

and D is the DMPK operator. Given the probabihty distribution Pl the conduc- 
tance of the wire can be obtained by means of the Landauer formula [15]. It is im- 
portant to stress that the construction of the DMPK equation requires an isotropy 
ansatz among the channels. Physically this implies that the finite time required for 
the electron to be scattered in the transversal direction is assumed infinitesimally 
small, and this assumption fits only for a quasi- ID wire. 

The Transfer matrix ensembles which are at the basis of the DMPK equation 
are very interesting from a mathematical point of view since they can be mapped on 
suitable symmetric space of negative curvature (see [2] and [3] for a detailed dis- 
cussion). Similar representations in terms of symmetric spaces also exist for other 
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Random Matrix Theories, but they usually involve symmetric spaces of zero or pos- 
itive curvature (see [3]). The negative curvature of the space is a peculiar feature 
of the transfer matrix ensembles and has relevant consequences on the quantum 
wires applications of these ensembles. For instance it can be shown that the eigen- 
parameters {Aj} appearing in formula (2) inside A exactly correspond to the radial 
coordinates of the symmetric space while the "radius" of the space can be mapped 
(with a suitable normalization) onto the length of the wire. The fact that the radial 
coordinates in a negative curvature symmetric space flow to infinity as the radius 
increases then becomes in the quantum wire context the well known result that as 
the length of the wire increases the wire moves from a conducting to an insulating 
regime. We shall discuss this mapping in great detail in the next section. 

Table 1 contains the main symmetry classes associated to the parameter (3 and 
the corresponding symmetric space X associated to the eigenparameters {Aj}. 



Tabic 1: Symmetry classes for a quantum wire 



/3 


TRS 


SRS 


M 


X 


1 


Yes 


Yes 


Sp{2N, R) 


Sp{2N, R)/U(N) 


2 


No 


Yes 


SU{N, N) 


SU{N, N)/SU{N) ® SU{N) ® U{1) 


4 


Yes 


No 


SO*{4:N) 


S0*{AN)/U{2N) 



3 Calogero-Sutherland models: mapping of the 
DMPK equation on a symmetric space and so- 
lutions 

In order to address the symmetric space description of quantum wires we need 
an intermediate ingredient, i.e. the well known Calogero-Sutherland models (for a 
review see [14]). These models describe N interacting particles in one dimension, 
and are characterized by an Hamiltonian that can be mapped, for certain values of 
the coupling constants, into the radial part of the Laplace-Beltrami operator [B in 
the following) on a suitable symmetric space X. This operator is the generalization 
of the familiar Laplace operator describing free diffusion in the Cartesian space and 
it describes free diffusion on a generic manifold X (in our case X is one of the non- 
compact symmetric spaces X of Table 1, characterized by a root lattice structure 
of type Cjv). 
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The action of the Laplace-Beltrami operator can be formulated in terms of the 
following eigenvalue equation: 

B<^k{x) = k'<^k{x) (4) 

the eigenfunctions ^ki^) are known in the literature as "zonal spherical functions" 
(ZSF) (for a review see [3]). 

The connection between the C-S Hamiltonian and the symmetric space (in par- 
ticular the root lattice associated to the symmetric space) and the mapping of the 
C-S Hamiltonian into the radial part of the Laplace-Beltrami operator is the key 
point which determines the integrability of the model. 

The general form of the Calogero-Sutherland Hamiltonian is: 

1 " 

d 

a; = (xi, ...,a;Ar), = -i — , = (x, a) (5) 



where A"^ is the set of positive roots a associated to the space X, (x, a) denotes 
the scalar product, v{xa) is the interaction potential and is a coupling constant 
which depends on the roots as: 

2 ma(ma - 2)|q;2| 
9a = ^ (6) 

where is the root multiplicity in the space X and |a^| its length; relation (6) is a 
necessary constraint for the model to be solvable. The particular C-S Hamiltonian 
we are interested in contains a potential of the form: 

v{Xa) = — (7) 

smh Xa 

The relevant point for the present analysis is that if we choose, as in [4] : 

Aj = sinh^ Xi (8) 

the DMPK operator can be mapped exactly onto a C-S Hamiltonian with a potential 
of the type (7) and the symmetries of one of the symmetric spaces listed in Table 1. 

These spaces are characterized by root lattices (denoted as Cjv) composed both 
by ordinary roots (whose multiplicity nio = 1, 2 or 4 can be identified with the /3 
parameter which encodes the symmetry properties of the matrix model) and long 
roots with multiplicity mi = 1 which arc responsible for the peculiar properties of 
the eigenvalues near the boundary A ~ (see [2, 3] for further details). Given the 
set {cj} of linearly independent versors associated to the space, ordinary roots are 
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commonly written as a = ie^ ± Ci {j ^ i — 1, N) and long roots correspond to 
a^±2ej {j = l,N). 

These root lattices are characterized by a Z2 symmetry, i.e. they are invariant 
under reflection. This means that the origin plays a distinguished role. Indeed, 
most of the results we shall discuss in the following can be more easily understood 
by introducing an additional (fictitious) eigenvalue which is kept fixed in the origin 
and interacts (with a standard repulsive interaction mediated by the long roots of 
the lattice) with all the other eigenvalues. 



Table 2: Symmetric spaces associated to the transfer matrix and their root multi- 
plicities 



p 


X 


ruo 


m 


rus 


1 


Sp{2N, R)/U(N) 


1 


1 





2 


SU{N, N)/SU{N) SU{N) ® C/(l) 


2 


1 





4 


SO*{4:N)/U{2N) 


4 


1 






The complete serie of steps involved in this mapping can be found in reference [2, 
3] . In particular it turns out that the radial part of the Laplace-Bcltrami operator 
B on the symmetric space X is related to the DMPK evolution operator D in 
equation (3) by: 

D = ^e{x)Br\x) (9) 
where the function ^{x) is given by: 

^{x) = I sinh^ Xj — sinh^ 1 2 | sinh 2xi \ 2 . (10) 

i<j i 

and the operator B depends only on the radial coordinates {xi} of the space ac- 
cording to the relation: 

As a consequence, if ^k{x), x = {xi, ■ ■ ■ ,Xj\f}, k = {ki,---,kj\f} is an eigen- 
function of B with eigenvalue /c^, then ^(a;)^$fe(a;) will be an eigenfunction of the 
DMPK operator with eigenvalue k'^/{2'y). 

Once the eigenfunctions of the DMPK operators are known it is rather straight- 
forward to construct the Green function and use it to solve the DMPK equation. 
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The solution turns out to be rather involved, but it drastically simplifies in the two 
interesting limits of insulating (small values of k) and metallic (large values of k) 
regimes. It is important to stress at this point that even if an explicit form is known 
only for the special case /5 = 2, a general and powerful asymptotic expansion for 
large values of x exists for all the possible symmetric spaces [2, 3]. This asymptotic 
expansion turns out to be enough for the construction of the two (insulating and 
conducting) limiting solutions. Its explicit form is: 



$fe(a;)~-^(Ec(r^K^^'=''M > (12) 



where rk is the vector obtained acting with r E W on k, W is the Weyl group 
associated to the root system of the symmetric space and {k, x) denotes the scalar 
product. 

All the informations related to the underlying symmetric space are encoded in 
the function c{k) which is given by: 

c{k) = n (13) 

a6A+ 

with 

c (k) = (14) 

where F denotes the Euler gamma function, a is a generic root belonging to the 
root lattice which defines the symmetric space, denotes its multiplicity and the 
product is restricted to the sublattice A+ of positive roots only. 



4 "Clustered" solutions of the DMPK equation 

In the previous section the integrability of the DMPK equation has been stressed 
in order to write explicitly, at least in the asymptotics, its solutions. In several 
physical applications (see next section) it would be interesting to study solutions 
of the DMPK equation in which the symmetry among the N eigenvalues is broken. 
The aim of this paper is to show how to exploit again the integrable nature of the 
DMPK equation in order to go beyond the isotropic solutions obtained in [2]. The 
key point is that, given the integrability of the DMPK equation, its solution can 
be written even if non trivial initial conditions are assumed; this observation allows 
to construct non isotropic solutions starting from an equation which is isotropic 
by construction. The simplest way to obtain this is to impose that some of the 
eigenvalues form a "cluster" (i.e. the distances among them remain finite while the 
distances with respect to the other ones go to infinity), or more generally a set of 
independent clusters; the peculiar form of the DMPK equation ensures that such 
solutions exist and survive in the asymptotic limit. These solution obviously require 
suitably chosen initial conditions, which are degenerate in the whole phase space. 
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The main goal of this paper is to show that also these clustered solutions admit 
an asymptotic expansion similar to that of cq.(12). Given this expansion, one can 
then obtain the probability density P{{xn}, s) in presence of these clusters in a way 
analogous to the isotropic case. We shall consider below some possible applications 
of this result. 

Let us see these solutions in more detail. Let us assume the cluster to be 
composed by the first N' < N eigenvalues. This means: 

\xi-Xj\<oo, iJ = l,...,N' {i^j)- (15) 

\xi-x,\^^, i = l,...,N- j = N'+l,N (ij^j) 

In the symmetric space framework we can identify the cluster by selecting a subset 
of the root system associated to the space. Let 11 be the system of simple roots 
associated to X, and 11' be a subsystem of simple roots which satisfies the inequality: 

n' = {a e n/ lim Xa < oo} (16) 

where 

There are at this point two possibilities. Since 11 is a Cat type lattice then 11' 
can be again of type Cn (in this case it must also contain the long root and the 
ordinary roots must be chosen so as to keep the Z2 symmetry of the lattice) or it 
can be of type Aj^. In both cases from the ordinary roots of 11' one can construct 
the differences Xa = Xi— Xj+i which correspond to the nearest neighbours distances 
between some of the eigenvalues. From the definition of 11' these distances must 
remain finite in the asymptotic limit so that 11' defines a cluster (if it is connected) 
or a set of clusters otherwise. If the cluster is of type there is no other costraint 
and the cluster can in principle fiow to an infinite distance from the origin (while 
keeping a finite distance among the eigenvalues inside the cluster). If the cluster is 
of type Cjv on the contrary the eigenvalues are bounded (by the lattice structure of 
n' itself) to stay within a finite distance from the origin. In the following we shall 
denote with x the set of radial coordinates outside the cluster and with x' the ones 
inside the cluster. 

The asymptotic expansion of the ZSF's in presence of such a cluster was ob- 
tained a few years ago by Olshanetsky in [16]. It turns out to be a rather natural 
generahzation of the Harish- Chandra result of eq.(12): 

*fe(x)~ ^ c,(rA;)e'('-''='*)*(,fey(x') (17) 

reW/W 

where "^ki^) is given by the product ^{x)^k{x). In this formula (rk)' denotes the 
projection of the vector rk on the sublattice 11' and rk its complement, ^{x) is given 
by eq.(lO) and 

c,{k)= n (18) 

aeA+/A'+ 
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where A'+ is the set of positive roots associated to the cluster, and the W which 
appears in the coset W/W is the Weyl group associated to the cluster. This 
expression is apparently simple but it is highly non trivial. Notice for instance 
that the symmetrization with respect to the Weyl coset W/W acts not only on the 
part containing the coordinates x but also on the momenta of the zonal spherical 
function describing the cluster coordinates ^(rk)' i^') ■ This means that the particles 
inside the cluster do not move independently in a section of the whole space but 
they "feel" the presence of the other particles and are subject to the symmetry 
group of the remaining space. It is interesting to notice that the above construction 
could be formulated also in the framework of the original Transfer Matrix ensemble 
(i.e. before diagonalizing the transfer matrix). The clustered solutions correspond 
in this framework to peculiar boundary limits of the original symmetric spaces 
(which are known in the mathematical literature as Martin boundaries [17, 18]). 
The precise characterization of these boundaries is outside the scope of the present 
paper. We plan to address this issue in a future publication. 

Let us see two examples which may hopefully clarify the issue: 

Example 1: a Cat type cluster composed by the first (A^ — 1) eigenvalues plus 
an isolated eigenvalue which flows to infinity, i.e. (see eq.(16)) 

\Xi-Xj\<00, i, j = 1, . . ., - 1 (iy^j). (19) 
\xi — xn\ oo, i — 1, . . . , N — 1 

Then the asymptotic form of the ZSF, and of the related product function ^k{x) 
is: 

N 

*fc(^)~Ec.(fc,)*.r(^')e^'^"" (20) 

where \E'^-^ (x') is associated to the ZSF of the Ctv-i symmetric space of the cluster, 

kj = (/ci, ■ ■ ■ kj^i, kj+i, ■ ■ ■ kj^f) denotes the collection of (N — 1) momenta in which 
kj is omitted and Cz{k.j) is given in this case by the product of the Ca{k) functions 
over all the (positive) roots of type a = iiej ±ei (V/ j), a = 2ej (j fixed). From 
this expression it is straightforward to construct iteratively the ZSF in which two 
or more eigenvalues flow to inflnity. 

Example 2: a Cat type cluster composed by the first two eigenvalues: xi and 
X2, while the other (A^ — 2) flow to inflnity. 

^k(x) ~ ^ c,(A;)*,(i),,(2)(xi,X2)e'S^=3*^^(")^'- (21) 
sew 

where W = W/W is the coset of Weyl groups associated to the whole set of eigen- 
values and to the cluster respectively, the subscript of the two-particle wavefunction 
^s(i),s(2)(a^i, 3^2) identifies the two momenta to which this function is associated and 
Cz{k) in this case is the product of the Cq,(A;)'s over all the positive roots of type 
2es(m), 2es(„) (with m,n > 2) and ±es(^) ± Cg^n) (with the only exclusion of the 
combination {m,n) — (1,2) or (2, 1) ). 



10 



5 Applications 



Among the various possible uses of these new solutions we see in particular four 
interesting applications. 

1] We can use them to model systems in which the number of open channels is 
reduced by the structure of the wire itself (see the wide-narrow-wide geometry of 
[19])^. In this case, one could consider a configuration formed by a Cjv type cluster 
(bounded to the origin) made of N' eigenvalues and let the remaining N — N' 
eigenvalues flow to infinity (see example 1 above). One can choose N' and s (the 
reduced length of the wire) so as to keep the cluster in the metaUic regime, while 
the other eigenvalues are in the insulating one and do not contribute to the wire 
conductance. 

2] The same configuration discussed above could be used to provide a simple 
but effective way to improve our description of real quantum wires. While in the 
previous example A'^' was fixed (being the number of the open channels in the 
narrow part of the wire) one could easily generalize the example keeping N' as an 
additional degree of freedom which can be varied so as to take into account the 
effect of external parameters like the amount of disorder in the wire. The rationale 
behind this proposal is the well known idea [20] that in a generic conductor only a 
fraction (which depends on the disorder) of its channels is open and the isotropy 
anzatz should be a very good approximation for them, while the non-trivial (i.e. 
non one dimensional) behaviour of the conductance is due to the variation in the 
number of open channels. 

3] The previous picture could be refined choosing configurations with more than 
one cluster, each cluster being characterized by a different distance at which the 
crossover between metallic and insulating regime occurs (fixed by the only free pa- 
rameter which we have in these generalized solutions, that is the number of eigen- 
values of each cluster). This multi-cluster structure could be used in order to evade 
the isotropy ansatz (this was indeed the main reason for the present investigation). 
In order to obtain the multi-cluster configuration it is sufficient to choose a not 
connected subset of roots 11' which satisfies the condition given in eq. (16). The 
roots belonging to 11' correspond to a set of coordinates x' which are grouped into 
different clusters; the number of clusters Nc and their width are uniquely deter- 
mined by the choice of the not connected root lattice 11'. The asymptotic expansion 
given in the general form of eq. (17) can then be applied in order to find the proper 
solutions of this multi-cluster configuration. 

The isotropic symmetric space characterized by N radial coordinates in this 
case is mapped into an anisotropic model in which the Nc channels correspond 
to Nc clusters (each cluster containing a different number N^ of coordinates); this 
reconfiguration should ensure a different probability for the electron to be scattered 
among the different channels. It is important to stress once again that in our 

^We thank C.W.Beenakker for suggesting us this possibihty. 
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proposal this goal is reached by properly choosing the initial conditions, while the 
DMPK equation (with all its remarkable properties) is kept unchanged. This means, 
as a side remark, that the DMPK equation still depends on only one parameter i.e. 
the ratio N/s. The Nc parameters which appear in the special solutions in which 
we are interested (see the general solution and the examples discussed in sect. 4) 
come from the initial conditions and simply consist in the list of the sizes of 
the A''^. chistcrs. In this type of solutions the channels will encounter the metal- 
insulator transition at different lengths of the wire; this feature suggests that the 
N^/s parameters can be related to the different mean free paths associated to the 
channels. This picture reminds the "non equivalent channels model" proposed by 
Mello and Tomsovic in [10] but differs from it since our proposal (as mentioned 
above) does not require to modify the DMPK equation. 

4] An interesting independent application is related to the recent attempt [8, 9] 
to describe the crossover between the metallic and insulating regime in quasi ID 
systems by separating out the first [8] or the first two [9] eigenvalues and considering 
the rest as a continuum in the solution of the DMPK equation. This approximation 
reminds the clustered solutions that we have discussed in this paper (see in par- 
ticular example 2). The fact that these clustered configurations are indeed exact 
solutions of the DMPK equation (albeit with non trivial initial conditions) may ex- 
plain the remarkable stability of the saddle point solution noticed in [9] and could 
offer an independent justification for the approach. 

Acknowledgements This work was partially supported by the European Com- 
mission TMR programme programme HPRN-CT-2002-00325 (EUCLID) . 
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